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ABSTRACT. We continue our exercises with the universal R-matrix based on the Kho- 
roshkin and Tolstoy formula. Here we present our results for the case of the twisted 

(2) 

af fine Kac-Moody Lie algebra of type A 2 . Our interest in this case is inspired by the 

(2) 

fact that the Tzitzeica equation is associated with A 2 in a similar way as the sine- 
Gordon equation is related to A| . The fundamental spin-chain Hamiltonian is con- 
structed systematically as the logarithmic derivative of the transfer matrix. L-operators 
| of two types are obtained by using ^-deformed oscillators. 
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J-l i 1. INTRODUCTION 

■ 

The Izergin-Korepin model [14] was introduced as a quantum integrable model 
related to a classical integrable system described by the nonlinear differential equation 

= -m 2 [exp(-2F) -exp(F)] 

for a function F of two independent variables. Although it was introduced for the 
first time by Tzitzeica within the framework of differential geometry, this equation is 
mostly known as the Dodd-Bullough-Mikhailov or Jiber-Mikhailov-Shabat model, 
since it was investigated later in the context of the theory of solitons and general as- 
pects of classical integrability. It is also known that this equation is a particular case 
of the Toda equations associated with twisted loop groups [24, 25]. Here the affine 

group of type A 2 plays the role of the underlying symmetry group, which for the 

sine-Gordon equation takes the simplest affine group of type A\ . 
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To a certain extent, the model under consideration has features reminiscent of both, 
the so-called 5(3 and sii, cases: on one hand it is related to a three-dimensional matrix 
representation, while, on the other hand, the infinite-dimensional part of the whole 
representation should correspond to only one scalar field entering the integrable equa- 
tion written above, and therefore we expect that a single copy of the ^-deformed oscil- 
lator algebra will be sufficient to describe L- and Q-operators. 

We work in the spirit and use notations of [5], continuing our exercises with the 
universal R-matrix based on the remarkable formula presented by Khoroshkin and 
Tolstoy [18, 19, 30]. After reproducing the ^-matrix of the Izergin-Korepin model, we 
construct two types of the L-operators using a ^-deformed oscillator algebra. In the 
case under consideration these L-operators turn out to be related by a similarity trans- 
formation. A similar relation holds for the K-matrix as well. We also discuss various 
(anti)automorphisms generating new L-operators from given ones. The L-operators 
can then be used for the subsequent construction of Q-operators encoding the infor- 
mation about the spectrum of the quantum integrable system under consideration. 

The problem of the calculation of the eigenvalues of the transfer-matrix for the 
Izergin-Korepin model was studied in [26]. An algebraic Bethe ansatz solution of 
this model was constructed in [27, 29]. This consideration was extended in [23, 22] to 
the case of open boundary conditions. 



2. Generators and roots 

t A = («i/)i,;=o,i be the 

form 



(2) 

Let A = (fljy)j ,;=o,i be the generalized Cartan matrix of type A { 2 ' having the explicit 



A 



2 -1 
-4 2 



The matrix A is symmetrizable and we have d ? a^ = diau for do = 2 and d\ = 1/2. We 
denote 



4 -2 



A = ( a ij) = ( d r a ij) =■ 1 _ 2 1 

Before proceeding to the universal R-matrix, let us first describe the structure of 
the Kac-Moody algebra g'{A), enveloping algebra U(g'(A)) and its quantum defor- 
mation Uf,(g'(A)). In accordance with the general construction, see, for example, the 
books [17, 10], the twisted affine Lie algebra g'(A) is generated by the elements hi, e{, 
fi, i = 0, 1, with the defining relations 

[hi,hj\=0, (2.1) 
[hi,ej\ = ciiff, [Kfj] = -<hjfy ( 2 - 2 ) 
[e u fj] = Sijhi, (2.3) 
satisfied for all i and /, and the Serre relations 

(adci) 1 "^*;) = 0, (ad/,) 1 -^) = 0, (2.4) 

satisfied for all distinct i and /. 

We denote the linear span of the generators hi by t)'(A) and its dual space by f)'*(A). 
The vector space t)'(A) is the Cartan subalgebra of g'(A), and we have the decompo- 
sition 

tf(A) = f)'(A)e © g'(A) Jf 

-yeA(A) 

where for any 7 £ t)'*(A) we denote 

g'(A) 7 = {x e g'{A) | [h,x] = 7(x)xforanyfr e f)'(A)} 
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and 

A(A) = {7 G It* (A) I 7 7^ 0, g'(A) 7 ^ {0}}. 

The elements of A (A) are roots, and g'(A) 7 is the root space of 7 whose nonzero 
elements are root vectors. 

The generators are evidently root vectors. We denote the corresponding roots by 
0L{. These are simple roots. Any other root is a linear combination of simple roots 
with integer coefficients all of which are either non-negative or non-positive. In the 
former case we have a positive root and in the latter a negative one. In particular, 
the generators ft are root vectors corresponding to the negative roots — aj. One has 
A (A) = A+(A) U A_ (A), where A+(A) is the set of positive roots and A_(A) = 
-A + (A). 

The symmetrized Cartan matrix A s determines a symmetric bilinear form on f)'*(A) 
by the equality 

{cci,ccj) = afj. 

Explicitly we have 

(ct ,ci ) = 4, Oo,ai) = -2, (ai,ao) = -2, = 1- 

•1 

It is convenient to denote 

5 = OCq + 10L\, Oi = Di\, 

so that 

(6,6) = (6,0c) = (a, 6) = 0, (tx,tx) = 1. 
It can be shown that the set of the positive roots is 

A+(A) = {a + md \ m e Z> } U {la + (2m + 1)6 \ me Z> } 

U {md I m G Z >0 } U {6 - lot + 2m6 \ m G Z> } U {6 - a + md \ m G Z> }, 

see, for example, the book [10]. All root spaces corresponding to positive and negative 
roots are one-dimensional. Choosing a root vector for each root and adding the Cartan 
generators hj we obtain a Cartan basis of g'(A). 

The enveloping algebra U(g'(A)) is defined as the unital associative algebra with 
generators hj, e,, /, and with the same relations (2.1)-(2.4) as g'( A). Here we can rewrite 
the Serre relations (2.4) as 

E 7 (-l)*^ 7") (^-^S (*i) k = 0/ (2-5) 

The Lie algebra g'(A) can be naturally considered as a subspace of LT(g / (A)), and a 
Cartan basis of fl'(A) as any of its bases generates a Poincare-Birkhoff-Witt basis of 
U(Q'(A)). 

Let Si be an indeterminate and q = exp h. In accordance with the general definition, 
see, for example, the book [11], the quantum deformation of LT(g / (A)), the quantum 
group 14 (g' (A)), is the topological C [[h]] -algebra generated by six elements hj, e,, f ir 
i = 0,1, with the defining relations 

[hi,hj]=0, (2.7) 



^In fact, 8 is the minimal positive imaginary root [17, 10]. 
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[hi, ej\ 



[hufj] = -<Hjfj> 



qdihi _ q-dihi 



-di 



(2.8) 
(2.9) 



satisfied for all i and /, and the Serre relations 

l-fl;; 



E (-i) fc 

k=0 

1-fljy 

E (-i) fc 

k=0 



{eif-^e^etf 



0, 



J f 



(0 



: fi(fi) k = o, 



(2.10) 



(2.11) 



satisfied for all distinct i and /. Here we have introduced the ^-binomial coefficients 



[n - m] q l[m] q V 



[nUn -!]«••• [l] q , [? 



It is clear that the relations (2.7)-(2.9) are deformations of the relations (2.1)-(2.3), 
and the quantum Serre relations (2.10), (2.11) are deformations of the Serre relations 
(2.5), (2.6). Khoroshkin and Tolstoy use a slightly different definition of a quantum 
group [30, 19]. We come to this definition using the rescaling of the generators 



hi — > d- 1 h j , 



q-q- 



q d i — q~ 



1/2 



fi 



q-q 



q d i — q 



1/2 



fi- 



After this, the defining relations (2.7)-(2.9) take the form 

[hi,hj] = 0, 
[hi,ej\ = affy [h u fj] = -afy, 



[eufjl 



q h < - q-* 1 ' 

ii —1 ' 



(2.12) 
(2.13) 

(2.14) 



while the Serre relations (2.10) and (2.11) preserve their form. Below we work in terms 
of the rescaled generators. Note that the element c = Hq + 2h\ belongs to the center 
of Uf,(g'(A)). It is convenient to assume that the definition of Ufi(g'(A)) includes an 
additional relation 

c = h + 7h\ = 0. 

This allows us to use for the case of the quantum groups LZ^(g'(A)) the Khoroshkin- 
Tolstoy formula for the universal R-matrix valid for quantum groups Uf l (g(A)) just 
putting there c = 0. 

The quantum group 11% (g' ( A) ) is a topological Hopf C [ [h] ] -algebra with correspond- 
ingly defined comultiplication, counit and antipode, see, for example, the book [11]. 

Associate now with each root from the infinite root system A (A) a corresponding 
root vector. We denote the root vector corresponding to a positive root 7 by e 7 , and 
the root vector corresponding to a negative root — 7 by f T Let us start with the simple 
positive roots and put 



-a 



e5-2a = CO, 



-a 1 
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The higher root vectors corresponding to positive roots are defined recursively by the 
relations 

es-a = ([2]q) 2 [eoc,es-2oc]q, e's = l e <x,es-a]q, (2.15) 

Cci+mS = ([3]gl/2)~ 1 [e a +( m -l)<5/4]<7' e 5-ot+m6 = ([3] (? i/2)" 1 [4/^-a+(m-l)<5)]'?' ( 2 -16) 

e 2ix+(2m+l)S — ([2]q) 1/2 [ e a+mS> e a+(m+l)s\q' (2.17) 

[ e S-a+{m+l)S> e S-a+ms]q> (2.18) 

e 'mS ~ [ e a+{m-\)S' e 5-a}q> {2.19) 

where the ^-deformed commutator is defined as 

for any two roots 7 and 7' from the root system A+(A). These relations allow us 
to construct the root vectors corresponding to all the roots from the set A + (A). In 
fact, the expression for the universal .R-matrix given by Khoroshkin and Tolstoy [30] 
contains the root vectors e m g related to the root vectors e' m& by the equality 

(q - q~ l )e S {x) = log[l + (q - q- v )e' 6 {x)}, 

where 

e 's( x ) = E e 'mS x ~ m > e s(x) = E e mSX~ m - 
m>0 m>0 

To complete the set of root vectors we need to construct root vectors for negative 
roots from the system A_ ( A). We do it using the rule 

fry = CO(ey) = e-sy 

for any 7 £ A+(A). Here to is the Cartan anti-involution defined on the generators 
by the relations 

coQu) = h, co(ei) = f if to(fi) = e { (2.20) 

supplied with the rule co(h) = —h implying that co(q) = q~ x . Note that we use such 
normalization of the root vectors that for any 7 7^ mS we have 

where h 7 = YU mfc if 7 = m^. It appears that, as in the case of U(q'(A)), the root 
vectors corresponding to all roots from A (A) together with the Cartan generators hi 
generate a Poincare-Birkhoff-Witt basis of Uf,(g'(A)). 

3. Universal K-matrix 

We consider the case of a quasi-triangular Hopf algebra, see, for example, [11]. The 
corresponding universal K-matrix 1Z satisfies the Yang-Baxter equation 

n 12 n 13 n 23 = n 23 n 13 n t2 . (3.1) 

and serves as basic object for the construction and investigation of integrable models. 
In particular, let for any £ £ C x a representation <p^ of Ufc(g'(A)) be given. Then the 
parameter dependent R-matrix 

R(CilC2) = ni®n 2 (tt) 

We use primed notation for the root vectors corresponding to the roots mS because we redefine them 
below. 
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satisfies the Yang-Baxter equation of the form 

^12(^1^2) ^13(^1^3)^23(^2^3) = ^23(^2^3) Ri 3 (^i|^ 3 ) ^12(^1^2)- 

Usually, one chooses the representation in such a way that R(£i [£2) depends only on 

combination C1C2 • This allows one to introduce the R-matrix R(z) depending on a 
single parameter, so that 

R(Ci|&) = R(CiG" 1 )- 

Now the Yang-Baxter equation takes its original form 

MC12) Rl3(Cl3) ^M&O = ^23 (£23) Rl3(Cl3) Ru(Cu)- (3-2) 

Here and below we use the notation = . 

An explicit construction of the universal R-matrix for quantum groups was pro- 
posed by Khoroshkin and Tolstoy [30]. For the case under consideration it looks as 
follows [19]. 

First of all, we have to choose a normal ordering [21, 1] of roots from A+(A). In 
fact there are only two possibilities. We use the ordering where the roots a + m5 and 
2a + (2m + 1)5 come first as 

a, 2a + 5, cc + S, 2a + 35, a + 25, 2a + 55, a + 35, 2a + 75, ... , 

then the roots m5 come as 

5, 25, 35, 45, 55, 65, 75, ... , 
and finally the roots 5 — a + m5 and 5 — 2a + 2m5 come in the order 

... , 5 -2a + 65, 5-a + 25, 5 -2a + 45, 5-a + 5, 5 -2a + 25, 5 - a, 5 -2a. 
Gathering, we can write 

7 + k5 -< mS -< (5 - 7) 4- 15, 

where 7 = a, 2a. Another ordering is the reverse to this one. 

Now we have all the ingredients needed to construct the universal R-matrix. Ac- 
cording to Khoroshkin and Tolstoy, it has the form 

1Z = 7Z~(S 7Z~s TZ^s 1C. 
The first factor is the product of the ^-exponentials 

K lrm = exp^ ((q - q~ x ) e 7+mS <g) f 7+m sj , (3.3) 
where 7 = a, 2a, m G Z>q. Here we use the notation 

q, = q -ton) 
and understand the ^-exponential as the series 

x 2 x n 



exp (x) = 1 + x + 7^7 + . . . + 7^ + . . . 

q n -l 



where 



(n) q \ = (n) q (n - l) q ■ ■ ■ (2), (l) q , (n) q - ——. 

q 1 

The order of the factors in TZ^s coincides with the chosen normal order of the roots 
7 + m5. The second factor is 

K„ s = exp i(q - q~ v ) £ b' 1 e mS ® f mS \ , (3.4) 

\ m>0 J 
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where 

b m = Ml (q m - {-l) m + q~ m ) . (3.5) 

Note here that 

oc+mSr Cnslq = bn e oc+(m+n)5' 

as it used to be for the untwisted case [19]. The factor lZ y s is the product of the q- 
exponentials 

Hg-^m = ex P ?7 ((<? - I' 1 ) e (S-y)+mS® f(S- 7 )+ms) > ( 3 - 6 ) 

where 7 = a, 2a and m G Z>q. The order of the factors in IZ^s coincides with the 
chosen normal order of the roots (S — 7) + mS. Finally, for the factor /C we have the 
expression 

JC = exp (h h a (g) h K ) . 
4. Finite-dimensional representation 
Given £ 6 C x , define the three-dimensional representation of the quantum group 



lift (g' (A) by the relations 

n {e K ) = C Sl (£12 + E 23 ) , ^(e*-2«) = T° (t^) 172 E 31 , (4.1) 

<pz(hec) = E n - E 33 , ?d h S-2a) = ~2Eu + 2E 33 , (4.2) 

n(f*) = ^ ( £ 2i + £ 32) , = r So ([2]^) 1/2 E 13/ (4.3) 



where s ; , z = 0, 1 are some integers, i = 0, 1, and the 3x3 matrix units Ejy are defined 
as 

(Eijjmn = 8im$jn- 

The basic property following from this definition and used in what follows is given by 
the relation 

Ef/Efcz = SjkEil- 

Note that (pd^S-ia) + 2<p£(/*a) = 0. If we have an expression for <p^(a), where a G 
li?i (g' (A)), in order to obtain the expression for <p^(co(a)), we should simply take the 
transpose of <pr[a), also changing the deformation and spectral parameters as q to q~ l 
and £ to We denote this operation by D 3/ where the index corresponds to the rank 
of the matrix it acts on. We will use similar operations for matrices of different ranks. 
Now, using the recursive relations for the higher root vectors as given above and 



equations (4.1)-(4.3), we obtain 

n(*«+mf ) = 1~ m r +mS (("l) m E 12 + q- m E 23 ) , (4.4) 

nih+ms) = q m r Sl ~ ms ((-l) m E 21 + q m E 32 ) , (4.5) 

<Pfa-*+ms) = -I"™ ?~ Sl+mS ((-!) m+lE 21 + r m " 2 E 32 ) , (4.6) 

cpdfs-K+ms) = ~q m r (s " Sl) - ms ((-l) m+1 E 12 + q m+2 E 23 ) , (4.7) 

fd^ + (2m + l)s) = (-1)* + ([2],j) 1/2 E13, (4.8) 

^(/ 2a+(2m+ lM) = (-l)-^ + l r 2s a -(2m +1 )s ([2],)V2 E31/ (4 . 9) 

Wfo-2«+2m*) = ^C 8-281 ^ ([2] f ) 1/2 E3^ (4-10) 



'Here and in what follows, we use instead of the integers sq and si the integers s = sq + 2s 1 and si. 
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9dh-2oc + 2m5) = (-l) m q 3m r {s - 2Sl) - 2mS ([2] q ) 1/2 E 13r (4.11) 

where m = 0,1,2, .. .. For the primed positive imaginary root vectors we have the 
following expressions: 

n^'ms) = <r m+1 r s (i-i)™- 1 En + K-irr 1 - i~ m ~ l \ ^ + q- m ~ 2 e 33 ) , (4.12) 

(Pdf'ms) = H m ~ l C™ (("I)" 7 " 1 En + [(-l) m <? - H m+l ] E22 + q m+1 E33) • (4.13) 
This allows us to pass to unprimed quantities via the relation 

(q - q-^cp^esix)) = <^(log[l + (q - f 1 )^^)]), 

where 

e 's( x ) = E e 'ms x ~ m > e s(x) = E e mS x- m . 

m>0 m>0 

The corresponding expressions for the unprimed generators f m s can be found then by 
applying the Cartan anti-involution. After some calculations we find 

nM = l ^ mS ((" 1 ) m " 1 Ell + \{-l) m q- lm - q~ m ] E22 + q~ 3m E 33 ) , (4.14) 

n(fms) = ^C™ ((-l)™- 1 En + [(-l) m q 2m ~ q m ] E 22 + q 3m E33) , (4.15) 
where m runs over the set of all positive integers. 

5. R-MATRIX 

In this section, we construct the R-matrix corresponding to the representation cp^ de- 
fined in the previous section. As usual, the most cumbersome part of the calculations 
is about the factor <p^ <8) (pr^H^s)- We obtain the following diagonal matrix: 



^\ 2 )-\{q-^\ 2 



1 ' 1 n ® E n + 1 + ^^ 2 s E22 <8> E22 + E 33 E 33 



+ _^ 12 (En ® E 22 + E 22 <g) E 33 ) + * fg s (E 22 (8) En + E33 <8> E 22 ) 



Ell ® E33 + , : , H _£? E33 



1 - & 1 + ^ s i2 1 - <r 2 ei 2 1 + rtfc 



-11 



where 

^ Q q m -(- 1 ) mj rq~ mm ' 

A useful property of the transcendental function A is that 

\(qC) - A(-0 + \(q~ l Q = - log(l - C). (5.1) 

Further, using equations (4.4), (4.5) and (4.8), (4.9) we see that cp^ ® f^ 2 {7t Kr m) com- 
mute with cp^ <g> (p^ 2 (7^2oi,n) f° r any m and n, where TZ a ,m and 7£ 2a , M are factors in 1Z^$ 
corresponding to the roots a + md and 2a + (2n + 1)8, respectively. Similarly, using 
equations (4.6), (4.7) and (4.10), (4.11) we see that cp^ (8 f^ 2 (T^s-Zoc,m) commute with 
cp^ x ® cp£ 2 (7Zs- a ,n) f° r any m and n. Here TZs-2a,m and TZs-a r n are factors present in 
corresponding to the roots 8 — 2a. + 2m8 and 8 — a + n8, respectively. Therefore, 
we can rearrange the factors entering cp^ ® cpr 2 (1Z^s) i* 1 the definition of the universal 
R-matrix, so that the factors corresponding to the roots a + md will come first, and 
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then come the factors corresponding to the roots 2a + (2m + 1)5. Similarly, we can re- 
arrange the factors entering cp^ (g> q)£ 2 (lZ y #) in the definition of the universal R-matrix 
in such a way that the factors corresponding to the roots 5 — 2a. + 2m5 come first, and 
only then come the factors corresponding to the roots 5 — a + m5. The same useful re- 
arrangement turns out to be also valid for the matter of constructing the L-operators. 
Hence, in the case under consideration we can write 

where the factors entering this definition are given by the expressions 



m>0 

rv 



m>0 



and 



m>0 



m>0 



Using the expressions for (g) cp^lZa), cp^ <g> cp^Tlia) and ^ (g> cp^ z (TZs-2oc) f 
fCi ® ^ 2 (^-a) presented in Appendix A, we obtain by the respective matrix mul- 
tiplication 



I + % 1 (E12 E 21 + E 23 E 32 ) 



1 



12 



+ frfe «S En ® E32 + ji+fi^r ffi E23 » E 21 



+ 



12 J r 2s i 



£l 2 E 13 ® E 31 



and 

Hi 



yd 



1 - & 



E 12 + E 



32 



-23 J 



(q — q l )q 2 C\2 Sl 
E 2 i E23 — — — 7-=z E 32 



1 + ^1 



1 + 1^12 A 1 *J ^12 



-12 



+ 



s-2s a 
12^12 



'31 



-13- 



(i-C52)(i + r%) 

For the simplest part of the calculations, the last factor in the definition of the universal 
R-matrix, we find the following diagonal matrix 

<PZ\ ® ^ 2 (^) = 1 E n ® E n + E n ® E 22 + q~ x En ® E 33 

+ E 22 En + E 22 <g) E 22 + E 22 <g) E 33 

+ g" 1 E 33 E n + E 33 ® E 22 + q E 33 E 33 . 
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Now, multiplying all the factors in the given order, we finally obtain 

<Ph ® nz( n ) = e x Mu)-Hi-Hh) r(£ 12 ), 

where 

R(0 = q En ® E n + [E n <g> E 22 + c(£) E 12 (8 E 21 

+D 9 (fl(0 Ei 2 <g) E 21 ) + E 22 (8 E n ] 

cr(0 [E n ®E 33 + q b(C) E12 E 32 + c(£) E 13 E 31 + &(£) E 23 E 21 

+ Q 9 (q b(0 Ei 2 <8> E 32 + c(0 Ei 3 (8) E 31 + &(£) E 23 (g> E 21 ) 

E 22 (g> E 22 + E 33 E n ] 

+ p(0 [E 22 (g>E 33 + a(£) E 23 (g>E 32 

+n 9 (a(£) E 23 (g) E 32 ) + E 33 (g) E 22 ] + q E 33 (g) E 33 . 

Here we use the notations 



q+iq-l^q-q-i + q-^-q- 2 ?* 



d(0 



[l-C s )(l + q-K s 



^ l-q-H 5 ' [U (l-q-H s )(l + q-H s Y 

Choosing so = 1, Si = we recover the R-matrix presented in [19]. To compare with 
the original result of [14], one should apply certain similarity transformation and also 
adjust the respective parametrization, see, for example, [20, 23]. Here we have that 

nMQ) = d(0, Oi{p(0) = q- 2 p(0, fliMO) = f MO* 

and we take into account that we have 

9 (Eij <g) 

Emn) — Eyf (g) E nm . 

The explicit matrix form of the R-matrix for the considered representation is given in 
Figure 1. It is clear, in particular, that diagonalizing this matrix means an independent 
diagonalization of two its 2 x 2 sub-matrices with p and one 3x3 block with a. Here 

q 2 C — 1 q 3 C + 1 

we obtain three different eigenvalues {q, q^-j — — , q^-j — — }, having multiplicities 

q ~ £ q + £ 

{5, 3, 1}, respectively. 4 

The R-matrices corresponding to different values of s = So + 2si and s\ are related 
by a change of the spectral parameter and a gauge transformation [9, 5]. In the case 



"^Compare with the untwisted case associated with Ao , where one has only two different eigenval- 
ues {1, i_ rs } with multiplicities {6, 3}, respectively. 
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Figure 1. 

under consideration we have the following relation: 5 

* (s ' Sl) (£i2) = [G(Ci) ® G(fc)] [G(fi) ® Gfo)]" 1 , 

where 

/ C Sl \ 
G(0 =01 (5.2) 

V o o r si / 

and £12 denotes the ratio 

Note finally that the expression for the iv-matrix remarkably factorizes. Indeed, the 
expression for the universal K-matrix can be written as 

Observing that the relation 

holds in the case under consideration, we can thus write 

R = R+R R , 

where we have denoted the upper-triangular, diagonal and lower-triangular factors as 

_i 

R+(Cl2) = Hx ® n 2 ( n « n 2*) = I + fik ( £ 12 ® £ 21 + £ 23 ® E32) 



+ 



+ 



(g-g l)(g-l + (g + q l )g ^ 



Cl2 E 13 ® £31/ 



5 It is implied here that R(£) = R^ (Q, Rt 1 ^ (£ s ) = R( s '°) (£), and a similar convention will be used 
for the L-operators. 
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Ro(Cn) = e-^ti+^Zh)^ nz {n^ 5 K) =qE u ® E n 

H En <g> E 22 + ^ — w a" 1 F n £ 33 

r ^ , ! + ^12 r ^ , 1-^12 



+ 1 _ fl -27s £ 22 ® E n + E 22 ® E 22 + _' ^ E 22 ® E 33 

+ -2% wi T -3% \ <7 E33 (8 En + , E33 ® E 22 + 4 E 33 ® E 33 , 

(1 - 9 2 £f 2 )(l + 9 3 £f 2 ) 1 - g 2 ^ 



MC12) = ^ ® nii^Ks-toKs-cJC) = ^+ i_/s 2 ( £ 21 ® E12 + E32 ® £23) 

1 G12 

- ( v+~y e 21 « e 23 - 'v^C'" e 32 « El2 

1+^12 1 + 9 %12 

( ? - r i )( ? + r i +( t? - i -i) r i ^ 12 )r 12 2si „ 

+ (wylTT^S El3 ' 

respectively. Here we also have the Cartan anti-involution relations 

rj 9 (R + (0) = MC), «9(Ro(0) = r 2 «o(0/ 

implying that 

o 9 (r(0) = r 2 R(0- 

6. THE SPIN-CHAIN Hamiltonian 
The transfer-matrix of the system based on N sites is given by the relation 6 

T(C|£i, ...,&,) = tr (l?oi(C/ei)Ro2(C/&) • • • #on(£/£n)) , 

where the string of Rofc(£/ £fc) ac ts in Vq ® Vi ® . . . ® and the trace is taken over Vq 
in the given representation. The corresponding Hamiltonian is defined by the formula 

dT(Q 



H= T~\0 



where T{Q = T(£|l, . . . , 1), and allows one to obtain the ground-state energy of the 
system, see, for example, [31, 32]. 

Direct calculations show that the Hamiltonian in the case under consideration is a 
sum of four terms, 

H = #(12) + #(23) + #(31) + #(123)' 

where 

H _ 1 4- F^ +1 ) - n- 1 F^ f(' +1 ) - a F^ f(' +1 )" 

(12) — „_„-iZj ^12^21 ^^21^12 " ^11 ^22 H n 22 c ll 

1 H l=l 1 J 

rr _ 1 M') F( /+1 ) + F^ - F^ F^ +1 ) - a F^ F^ l+1 ^' 

n {23) — _ 1 23 ^32 32 23 1 n 22 n ?>3 1 n 33 n 22 

H H i=i 1 J 



6 Also a certain twist operator can be introduced in the definition of the transfer-matrix corresponding 
to specific boundary conditions. 
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H 



(31) 



[%/2 



N 



p(l) F (Z+1) p(Z) p(Z+l) 

„ n -l [ol L-i c 31 c 13 c 13 c 31 

q q p\q l=1 L 

-( r 1 +^-l))4 ) £j Z 1 +1) - 



(9 



r 2 to-D)ES4 +1) 



#(123) - 



[3/2], 



[3] 



E 



<7 Z=l 



1 F (0 F (/+l) 
32 ^12 



+^ 1/2 ^ ^ 1 E23 E21 ^) (*7 

in the standard notation, with s, fixed as So = 1,S\ = 0. 

Recall that for the untwisted case, i. e. the system associated with A K ^_ V we have the 
Hamiltonian also expressed in terms of the respective matrix units, 



1/2 _ --1/2 
n 



) [2] q E22 E 22 



(1) 



N 



H r , 



tE E E ij E 

--1 U,j=i 



In the best studied cases, n = 2 and n — 3, the Hamiltonian can be represented by 
means of the generators of the 1st fundamental representation of the corresponding 
finite-dimensional algebras, A\ and A<i, respectively, 



H 2 = 



1 



N r 



4 



4 1 / 4 q-q- 1 



N 



Ho 



+ 



q-q 



q + q- 



Z=l 



F (Z) F (Z+1) p(Z) p(Z+l) p(Z) p(Z+l) 
c a r a c a+(6 r a+ ^ -t- 

F (Z) p(Z+l) p(Z) p(Z+l) p(Z) p(Z+l) 



4K } 



(1) _ H (N+1)\ N_q±q_2_ 



3 q — q~ 



(2) 



Here we kept the boundary terms and the constant terms. 

Now, comparing these expressions with what we have obtained for the twisted A 2 ' 
case, we see that the Hamiltonian of the system under consideration is the sum of two 

pure -type systems' Hamiltonians H( 12 ) and H( 2 3), related to the indices (12) and 

(23), one more somewhat shifted A^ -type system's Hamiltonian H( 31 ), and an essen- 

(2) 

tially A 2 addition presented explicitly by the part Hm 2 ^\ . Here, we cannot express 

the Hamiltonian H in terms of matrices representing the algebra generators, which is 
in contrast with the untwisted cases. 
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7. L-OPERATORS IN ^-OSCILLATOR REPRESENTATION 

A useful object to investigate the properties of an integrable system is the corre- 
sponding Q-operator. According to the modern approach, it is constructed as the trace 
of some monodromy type operator constructed, in turn, from an L-operator. Here 
the L-operator is obtained by taking one of the factors of the universal R-matrix in 
an infinite-dimensional representation. Usually, it is some ^-oscillator representation 
[4,2,5]. 

In the case under consideration, 1Z is an element of the tensor product of the Borel 
subalgebras of the quantum group U%{q' (A)), and thus we have 

U G Uf l (b' + (A)) ® U h (b'_(A)) c U h (Q'(A)) <g> U Tl (g'(A)). 

The two Borel subalgebras here are unital associative algebras generated by the ele- 
ments hi, e\ and hj, /, respectively. To construct an L-operator it therefore suffices to 
have representations of the Borel subalgebras. 

7.1. Resolving the Serre relations. In this section we construct L-operators based on 
the ^-deformed oscillator algebra defined as an associative algebra Osc?j with genera- 
tors a, cr and D subject to the relations 

[D,a] = —a, = a + , 

~f~ o on -f- or) 

aa = 1 — q q , a a = 1 — q ■ 

The transformations 

a^Kaq^ D , a f ^K~ l q~^ D a\ D^D (7.1) 

form a two-parametric group of automorphisms of Oscf,. One can use these transfor- 
mations to obtain different L-operators, however, the trace used in the definition of 
Q-operators is invariant with respect to the action of this automorphism group [4, 7]. 

To construct L-operators we have to consider representations Xi, an d tfr °f (A)) 
and Ufi(b'_(A)), respectively, in addition to the representation cp^ already described 
in section 3. Slightly more abstractly, we will use homomorphisms Xl an d °f 
Ufi(b' + (A)) and U%(b'_(A)) to Osc^. It is easy to switch to representations, when nec- 
essary, using the well known representations of Osc^. 

By an L-operator of type L we denote an element of Osc^ (g> End(C 3 ) ~ Mat^Osc?,) 
defined as 

£(Cl2) = (ft). 

Here, we have assumed that the homomorphisms Xl an d <P£ are such that they result 
in an L-operator depending only on £12 = Cl/^2- It follows from the Yang-Baxter 
equation (3.1) that the L-operators of type L should satisfy the equation 

^23(^12) MCi) Mfc) = tufa) M&) MM- (7-2) 

As an L-operator of type L we define an element of End(C 3 ) ® Osc?j ~ Mat3(Osc?j) 

L(Ci 2 ) = n 1 ®^ 2 (ft) / 

^Instead of the naturally ^-deformed oscillators, defined by the relations [Af, b] = —b, [Af, = b + , 
bb* = [Af + l]q, b^b = [Af]q, we use slightly different objects used earlier in [6, 7, 8]. The latter are 
related with the former as D = A/", a = —(cj — q~ l )bq^ , a* = b f . 
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assuming again that the homomorphisms ipr and cp^ are such that they result in a 
dependence on the ratio Ci/Ci- Using the Yang-Baxter equation (3.1) we derive the 
following equation for the L-operators of type L: 

RniCu) iiaCCi) £23 (fc) = L 23 m £l3(Cl) Rl2(Cl2)- (7.3) 

In terms of the matrices R and R defined by means of the matrix of the permutation 
operator as 

A(0 = R(C)P, R(0 = PR(Q, (7.4) 

equations (7.2) and (7.3) take the forms 

A(C 12 )(L(Ci) H £(&)) = H t(Ci))A(Ci2) 

and 

£(Cu)(i.(Ci)BL(C2)) = (I(c 2 )HL(Ci))*(Ci2)/ 

respectively, where M means a generalization of the Kronecker product to the matrices 
with arbitrary algebra-valued entries [11, 5]. 

Remember that the Borel subalgebra Uf,(b' + (A)) is generated by the elements h ir e ; , 
while the dual Borel subalgebra Ufi(b'_(A)) is generated by hi, fj, i = 0,1. Here, the 
corresponding Serre relations (2.10), (2.11) are explicitly of the form 

4-2a e * ~ (<? 2 + 9~ 2 ) e <s-2« e« es-ioL + e* ej_ 2a = 0, 

5 rcl 4 , [ 4 L 1/2 [ 5 V /2 / 3 2 2 3n 

e a e,j-2a - P\^n e u e s-2a tot H (e« e<5-2 a e a - e a e ( 5_ 2a ej 

+ [5] q i/2e a es-2a A - es-iot e l = 0, 

and we have the same equations for fs- 2u and f u . Besides, the Cartan generators have 
to satisfy equations (2.13) and the condition hs-ia. + 27z a = 0. 

To fulfill the defining relations of (b' + (A) ), including also the corresponding Serre 
relations, we use the homomorphism x defined by the equations 

X(hs-2a) = 2D, xM = ~ D > 

X {e S - 2K ) = Mo « +2 r 2(l/+1) ° XM = Fiaq vD 

with free parameters fiQ, ]i\ and v. The ^-dependent homomorphism Xl can be defined 
by the same procedure used earlier for the homomorphism cp^. Here, changing the 
parameter jiq we always change the coefficient at £ s . Note that the parameters }Iq, 
v can be freely changed by the transformations (7.1). 

Analogously, to satisfy the defining relations of M^{b'_ (A)), meaning also the corre- 
sponding Serre relations, we use the homomorphism \p defined by the equations 

tp{h S -2«) = -2D, tp{h«) = D, 

Ms-*) = ?*o« + V 2(1/+1)D , </>(/«) = n*q vD > 

where jio, ]i\ and v are again free parameters. The ^-dependent homomorphism ipr can 
be defined by the same procedure used earlier for the homomorphism cp^. Changing 
the parameter fiQ we always change the coefficient at £ s , and the parameters jiQ, 
v can be freely changed by the transformations (7.1). In both cases, the L-operators 
corresponding to different values of y.Q, ]i\ and v are equivalent. 
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7.2. L-operators of type L. Thus, to construct L-operators of type L we use the homo- 
morphism xz, from Ut,(b' + (A)) to Oscf, defined by the following relations: 8 

Xd e S-2a) = 1 —=a n q- 2D ^, xd e «) = —^i"C Sl - ( 7 -6) 

(r-vJW, q - q 



This corresponds to the choice of the free parameters ]Iq = (q — 1) 1 ([2] £? ) 1 ^ 2 , }i\ = 
l/(q — q~ l ) and v = 0. 

The higher root vectors are defined according to the recursive relations (2.15)-(2.19). 
Applying the homomorphism Xl we subsequently obtain 

n m rms 

v (J \ _ 1 b r/-, n 2m-l\ /-, n 2m+l\ n 2D] _2(m-l)D 

K {6ms) - (q-l)(q-q-l) [(1 " 9 } " (1 ~ 9 " ] * ' 

that gives 

= - ^iFh) ( 1 -¥u Kqmq2mD )' 

where the quantities b m are given explicitly by (3.5). Taking also the expression for 
<pr(fms) from equation (4.15) we obtain for the image of the factor IZ^g the following 
expression: 

® VbV^S) - e I C 1 + 9 bl2<7 ) h H H , 4^s 2D ^22 

\ 1 + 9 bl2*7 

(i+, 2 a 2 )(i-, 3 a 2 ) 

The simplest part of the calculations is, as usual, given by the operator K.. In our 
case it is represented by a diagonal matrix of the form 

XZr ® <P&(£) = q~ D E n + E 22 + <? D E 33 . 
As we already did for the R-matrix in section 5, we can rearrange the factors entering 
Xl\ ® y&^^s) ^ n * ne definition of the universal R-matrix in such a way that the fac- 
tors corresponding to the roots a + m5 come first, and only then come the factors corre- 
sponding to the roots 2a + (2m + 1)5. Similarly, concerning the factor %g <g) cp^ijt^g), 
we rearrange the factors entering there in such a way that the ones corresponding 
to the roots 8 — 2a. + 2mS come first, and finally come those ones corresponding to 
5 — a + mS. 

Altogether, the relations derived in Appendix B allow us to write down expressions 
corresponding to the factors 1Z^$ = TZ K T^2a. and TZ^s = T^S-2uT^s-a.- We obtain 

x Cl ® n 2 (K^ s ) = i + ffi a (i + <? 2 C12 <? 2D ) _1 £ 21 



+ « (1 - ^ 3 a 2 ^ 2D )" 1 E32 + ^ $ ^ (1 + q* ^2 ^ 2D )" 1 E 31 



and 



Xlx ® n 2 (Kys) = I - q —^~ q 2 CrT 1 « + (1 + <? 4 & E12 



As above, also in what follows the more convenient combination of the integers s = Sq + 2s i and Sj 
will be used instead of the initial ones, Sq and Sj . 
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q - q~ ] 



ij-ij 4 s- Sl t ^5 r s n 2D\-\ F 



q- q s _ 2si +2 2D , „6 r s „2Dn-1 f 
q i 

Multiplying the above presented factors in the given order, TZ^$ 1Z~s TZ^s Kl> we finally 
obtain the L-operator. According to our terminology, the corresponding object is called 
an L-operator of type L, and therefore, we will use the respective notation. In matrix 
form we have 



m = e A (-^ s 



^aq~ D l-q? 1+1 a + q~ D 



We have the following relation between the L-operators with arbitrary values of the 
parameters and fixed ones: 



where the matrix G(Q was defined earlier by (5.2), while the mapping 77, £ € C x , is 
defined by the relations 

7f (a) =ap, 7^( fl+ ) = 7f(D) = D. 

It means, in particular, that the Q-operators based on L-operators with different val- 
ues of s and S\ are related by a change of the spectral parameters and a similarity 
transformation. 

Applying to L the automorphism a generated by the transformation 

a f ^aq~ D , a^—q~ D a f , D — > — D — 1, 
we obtain another L-operator of type L, 

a : L — > Lg-, 

where 

l,(o = e A (-<?n tua 

and we have used the notation 



K(0 



( q(q D +?q- D ) -(q + l)q C s ~ Sl a q' D (q + l)q £ s " 2s i a 2 q~ D \ 



V ^T^ 2Sl « t2 r D -q- l ^aU~ D q-\q- D + q 2 ? q D ' 



J 



Note also that an object similar to L' ff was used in [12] in attempts to describe so- 
called defects in affine Toda field theory. 

7.3. Type L and some useful relations. The L-operators of type L can be constructed 
by means of the homomorphisms cp^ and \pr- This requires to map the Borel subalgebra 
Uf l (b' + (A)) to the finite-dimensional matrix representation, while realizing the Borel 
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subalgebra Uf l (b'_(A)) in Osc^. Here we consider the following homomorphism for 
the generators f ao , f ai and h ao , h Ul : 

Ms-*) = — ^=« +2 ^ ^ 2D ^\ uu) = v^ a ^~ Sl 

(?-1)a/[2]« 



and 

M h s-2a) = -2D, ip c (h a ) = D, 

Using the Cartan anti-involution we obtain from the recursive relations (2.15)-(2.19) 
recursive relations for the higher root vectors spanning the whole Uf l (b'_(A)). Per- 
forming the whole procedure in the same way as in the preceding subsection, we ob- 
tain the following matrix form of the L-operator of type L: 

/ q D + C s q~ D pa ^Tl^ la2 r D \ 

i+l^-s la f q -D X _ qC s ^ aq -D 

yq_ + l C s-2s Xa f2 q -D _( q + 1)q p-* a i q -D + q 2 C s q D J 

Here, to bring the scalar factor to the given simple form, we have used the relation 



1(0 = e A ("^ 



M-qC) 



(7.7) 



(l + q 2 Q(l-q 3 Q 

following from the basic equation (5.1) for the function A. For the L-operators of type 
L we have the following relation: 

= G(d) 7& (^(Cu)) 

Now, comparing the explicit forms of the L-operators of the type L and L we see 
that they are related as 

where / is the 3x3 skew-diagonal unit matrix. 9 Besides, for the matrices R and R 
introduced by equation (7.4), we obtain the following relation: 

R=(J®J)R-(J®J). 

Also the inverse of the L-operator of type L at £ _1 should be an L-operator of type 
L. Explicitly we have 



where we have denoted 

/ q D -q 3 £ s q- D (q + l)q 3 ^ q~ D fl + -(q + l)q 3 C s ~ 2si q~ D a i2 \ 



< Sl a l + qH s -(q + l)q 3 ?- s ia f 

^q~ D a 2 -^q- D a q~ D - q 5 £ s q D 



l'(0 = 

And again, using equation (7.7) one can simplify the form of the scalar factor at V . 



9 This expression Jx J implies actually the superposition of two operations, where one has to take first 
the usual transposition of x with respect to its main diagonal, and then the transposition of the obtained 
matrix with respect to the skew diagonal. 
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Note finally that, if L(£) is an L-operator of type L, then t(L(£ -1 )) is an L-operator 
of type L, and vice versa, if L(£) is an L-operator of type L, then t(L(£ -1 )) is an L- 
operator of type L, where T means the anti-involution of the ^-oscillator algebra Osc?j 
defined by the relations 

r(a) = a + , T(fl f ) = a, t(D) = D. 

Besides, applying to L the automorphism a described in the preceding sub-section, 
one can obtain another L-operator of type L 

8. Concluding remarks 

We have constructed the R-matrix for the twisted Kac-Moody algebra of type A 2 , 
that is given by the universal R-matrix under the action of a twisted evaluation ho- 
momorphism. That latter allowed us to bring the dual Borel subalgebras to a finite- 
dimensional matrix representation. We have found an expression for the respective 
spin-chain Hamiltonian. We have also constructed L-operators realizing one of the 
Borel subalgebras in the ^-deformed oscillator algebra. 

We have seen, in particular, that the R-matrix and L-operators have two points of 
degeneracy in this twisted case in contrast to the earlier considered untwisted cases 
where we had only one such a point. 

Recall that one usually considers a useful decomposition of the R-matrix in the form 

R(z) =zR -z- 1 Rq\ 

where the non-degenerate matrix Ro does not depend on the spectral parameter z. 
Instead, taking, for simplicity So = 1, Si = 0, so that s = 1, we will have another 
relation, 

R(0 = A(0 [q l, \r l ri - r Q ) - 9 (q 1 ^ 1 n - r Q ))'_ , 

where the upper-triangular matrices tq and r\ do not depend on £ and the scalar pre- 
factor A is a rational function subject to the condition 

fli(A(0) = -q- 2 A{Q. 

Here we explicitly have 

r = En ®E n + q E n ® E 22 + E u ® E 33 

+ q (e 21 ® E n + (l-(q- q~ 1 ))E 22 ® E 22 + E 22 ® E 33 ) 
+ E 33 ®E n +q E 33 <S> E 22 + E 33 ® E 33 

-(.1- <T 1 ) (q~ 2 E 12 ® E 2 i - q E12 <S> E 32 - E 23 ® E 21 

+ (q + q~ l )q E i3 ® E 31 + q~ 2 E 23 ® E 32 ) , 

ri = q 2 En ®E n +q E n <8> E 22 + E n <8> E 33 

+ ^ (E 22 ® En + E 22 <8> E 22 + E 22 ® E 33 ) 
+ E 33 ® E u + q E 33 ® E 22 + q 2 E 33 ® E 33 

+ (q - q~ l ) (q E n ®E 21 + q e X2 ® e 32 + e 23 ® e 21 

+ (q-l) E B ®E 3 i + q E 23 ® E 32 ) 
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and 

Ail) = Z - 

We have also constructed L-operators of two types, L and L, in the ^-oscillator rep- 
resentation. They satisfy Yang-Baxter equations with the R-matrices R and R, respec- 
tively Here, the L-operators and R-matrices of different types are related by similarity 
transformations. 

The L-operators allow for the decomposition of the usual form £L + — £ -1 L_, where 
the constituents L + and L_ are presented by non-degenerate matrices. The latter dif- 
fers from the untwisted case, where one of such two parts, L + or L_, turned out to 
be degenerate [5]. However, to find also for the L-operators a decomposition simi- 
lar to what we have for the R-matrix, one needs to define an analogue of the Cartan 
anti-involution in the ^-oscillator algebra. 

An interesting and quite non-trivial question in a lattice model is about continuum 
analogues of any monodromy type matrices. In a reasonable classical limit one should 
anticipate that the corresponding objects, in terms of suitable phase space variables, 
would have certain fundamental Poisson brackets with the classical r-matrix. This 
problem, in a somewhat reversed form, was formulated in the pioneering papers on 
the quantum inverse scattering method, see, for example, [28, 15]. And also in our 
case it will be the matter of further investigation. 

Another special question concerns the quasi-classical limit. If we evaluate our L- 
operators naively at h — > 0, where q — > 1 and £ — > 1, we obtain a rather meaningless 
degenerate matrix. Instead, we could follow an idea exploited e. g. in [3] and first 
renormalize the L-operators by means of certain linear transformation. Let us consider 
L(£) of section 7.2 with fixed s = 2, si = and change there £ 2 by — £ 2 for convenience. 
Performing now the transformation 

where 

l q =( 1 ], % = JZ q J, 

we rewrite the transformed L-operator in terms of the naturally ^-deformed oscilla- 
tors 10 J\f, b and & + and use the parametrization q = e h and £ = e fa . Then, in the limit 
h — > we have 

/ z + Af + 1 2fc + 
t c (z) b 2 

V b 2 /2 b 

One can treat the obtained matrix L c (z) as a local quantum Lax operator for the quasi- 
classical limit of the Izergin-Korepin model. 

The objects constructed in this work will be used in future work to develop the 
method of functional relations including the transfer matrix and Baxter's Q-operators 
for the quantum integrable system considered here. Note that to obtain functional 
relations one can also use an approach based on the notion of fundamental modules 
over the Borel subalgebras of the quantum algebras [13, 16], allowing one to avoid 
explicit forms of the L-operators. 




10 See the footnote in section 7.1. 
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Appendix A. Calculating the R-matrix 

The calculation of the factors cp^ q>^ 1 {J^2a) an d <P£ a ® (p^ 2 (^s-2a) belongs obvi- 
ously to the simplest part of the work, since, due to the relation 

(E^ = 0, / / £ > 1, 

the infinite products in the corresponding expressions reduce to simple infinite geo- 
metric series. Indeed, we have 

rx 

Hi ® <PC 2 (^2 a ) = II \ l + fa ~ ni( e 2 a +(2m+l)^) ® (Pt 2 (f 2 a+(2m+l)s) 

m>0 

= /+(</- J] (P^(e 2c( +(2m+l)s) ® n 2 (/2«+(2m+l)j)/ 

m>0 

where J stands for the 9x9 unit matrix. Then we easily perform the summation in the 
expression above and obtain 

£S+2sj 

<PCi ® ffe^to) = I + (q ~ q' 1 ) [2] q E 13 ® E 31 . 

1 W2 

Similarly, for the factor lZs-2a we obtain 

<Ph <g> <p^ 2 {Hs-2a) = E[ (/ + fa _ <PCl( e <5-2a+2m<5) ® (P^ 2 {fs-2a+2ms) ) 

m>0 

= I+(q - q~ V ) Yj <PCii e 5-2a+2ms) ® 2a +2m<0 • 

m>0 

This allows us to write down the expression 

rs-2s! 

ni ® n 2 ( n s-2a) = i + (q- q~ l ) [2], % ® e b . 

The most complicated part of our calculations concerns the factors 7^ and 7^$—%. 
This is due to the fact that the images of the generators in the corresponding expo- 
nentials are now nilpotent of degree three. Hence, we need to take into account terms 
quadratic in such generators, 

rx rx I ^2 \ 

ni ® n 2 (^«) = n ex p,«(*«,»n) = n 1 + + t^vt ' 

while all higher order terms vanish. Here, for convenience, we have used the notation 

Xa,m = (q ~ q~ l ) (P^a+ms) ® <P ^ 2 {f oc+ms) , 

where ^^(e a+m j) and <p£ 2 (/ a+m< 5) are explicitly given above. The infinite product in 
the expression for cp^ <g> q^IZo) can be rewritten as the following infinite sum: 

m>0 q + i- m > m>0fc>0 
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Performing the corresponding summations, 



_1 

E Xcc,m = \ _\ s Cl2 ( E 12 ® E 21 + E 23 ® E 32 ) 

+ £& £32 + {£ E23 « En, 

7-1 > 



— T 2^ X «,™ - 1 721 *12 h 13 ® ^31/ 



we obtain 



1 



^ ® <p^ 2 (^ a ) = I + ^ ffi ( £ 12 ® E21 + E23 <8> E32) 



'12 



We have the same situation with the factor TZs-a- Here we have to carry out the 
expression 

® ni^-oc) = n ex P 9 «(y*,«) = n 1 + y*,« + 7^7 ' 

m>0 m>0 V y^JI*-/ 

where now we have denoted 

yot,m — (q - q~ V ) Cp^ies-a+ms) ® <P£ 2 {fs-oc+ms)> 

with <p£j (es-a+ms) and (f$- a+m s) given in the preceding section. Again, we rewrite 
the infinite product in the expression for cp^ ® (p^TZg-^) as an infinite sum, 

<PCi ® <Pz 2 {Ks-et) = 1 + E y<M» + 7X7 E ya,m "t" E E Voik+m+l }/cc,m- 

m>0 t m>0 m>0fc>0 

After some calculations, 

E y«,« = (? ( £ 2i ® £12 + e 32 ® e 23 ) 



m>0 



- 1 ,^ 12 E 21 ® E 23 - W ^- _Vs 2 £ 32 ® E 12 , 



-— T L y«,m - ; JTs E 31 <g> t 13/ 

V ^ 1 m >0 1 ^12 

E E y«,fc+m+l y«,m - _ n , .-Irs VI ?27T E31 ® El3 ' 

we obtain 
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cp Cl <8> (p C2 (Ks- a ) = I + {q ~ q _* 12 (£21 <8> E 12 + E 32 (g) E 23 ) 

r+SS 21 23 i+rt s 12 32 12 

(q- r %-i-(i- r V^ 12 )^ s - Si) 

(l-cgXi+rtfc) 

Appendix B. Calculating the L-operators 



With respect to the higher roots 2a + (2m + 1)5 we obtain the expression 

(,-l)[2]^.+P*+'>» , 2f;m+1 > D 

which, in addition to (4.9) for <p^(f2a+(2m+i)s)' leads to 

-v- m CV \ — T J- i(i~jO r s+2 Sl n 2 n 2D(-, j_ n 3 r 2s „4D\-1 tt 

Related with the roots £ — 2a + 2m5 we obtain the expression 

,8m 7"s — 2sj +2ms 



v ( a \ _ g b „t2 J2(2m-1)D 



(^-1)[2]? /2 



Afo <8> <P£ 2 (Kct) = I + E X *<m + — J] X l« + E E *«,fc+m+l/ 



J'? 

and so, using also equation (4.11) for ^Afs-ioc+ims)/ we obtain 

XCl n2 (^_ 2a ) = J + r i2 " 2si a +2 r 2 °(l + <? n C?2 El3- 

Further, for the factor TZ K we obtain the equation 

_g_ y 2 

m>0 " T m>0 m>Ok>0 

where now 

= (q - q~ V ) XCi( e oc+ms) ® Cp^ 2 (f a+m s)- 

Using the expressions 

/ \ H b ?mD 

X^x+ms) = ~ rr fli ? 

and (4.5) for fr(foc+ms)/ we obtain 

E = ffi fl K 1 + ^ & ^ 2D )" 1 £ 21 + (1 " <? 3 C?2 ^ 2D ) _1 E32], 

m>0 

g V v 2 — ^^12 „2 /-I , „3 r 2s n 4DN-l p 



m>Ofc>0 

This leads us to the expression 

XCl ® ffcfat) = I + f 12 a (1 + q 2 £{ 2 q 2D )~ l E 21 

+ eg « (1 - ^ 3 & ^r 1 e 32 + ^ & « 2 a + ? 3 & o~ a % 
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2 r s+2s l n 2 n 2D , „3 r 2s n 4D\-l(-> , „2 r s „2D\-1 p 
-q Cn a 1 \ l + C l Clll ) \ l + c l fel2<? J fc 31- 

For the factor 7^_ a we have the equation 

U\ ® = 1 + E y*/" + t-tt E yiU + E E j/a+m+i y«,m, 

m>0 t m>0 m>0fc>0 

where we have used the notation 

y« 7 m ~{q- q' 1 ) XCi( e S-a+ms) ® (P&ifs-a+ms)- 

Using the expressions 

xc\ e s-oi+ms) — ——, a q 

q 1 

and (4.7) for <p z {f s 

—a.+mS 

), we obtain 



E y*s» = -h^v? ST 81 fl K 1 + <r & 9 ) - <f a - r & <? 2D r £ 

m>0 9 1 

g V- 1 -.2 _ (^"^ X ) 2 r2(s-sj) +2 /-i , „11 r 2s _4D\-1 p 

^TT E/^""^! (^-1)2 ^ « (1 + 9 Cl29 ) Ei3, 



E E Vix,k+m+l Voi,m 
m>0k>0 

_ (q-q 1 ) 2 12 r 3s-2 Sl +2 „2D /-, , _11 r 2s n 4D\-l(-i , „6 r s _2D\-1 p 
~ (t?-l) 2 ^ 9 + tl2^ ) + Cl2<? J h 13- 

Summing up these terms, we obtain 

Xfc ® fa-*) = I - q 2 ^2 S1 ^ (1 + <? 4 & ^ D ) _1 El2 

q-q 4 s _ s + 5 _2D\-1 p 

+ q in a ^-q Cni ) h 23 
q 1 

<f {q-q' 1 ) 1 r 2{s- Sl ) +2 , 11 r 2s fl 4D\-l Fio 

, {q-q ) n 12 ^3s-2 Sl +2 „2D , „11 r 2s „4D\-l/-i , „6 r s „2D>>-1 p 
+ ^ ^2 9 6 12 fl 9 U+9 tl2 < ? J I 1 + 4 6 12 4 J h 13- 
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